home *** CD-ROM | disk | FTP | other *** search
/ AmigActive 26 / AACD 26.iso / AACD / Programming / ace_gpl_release / src / dp / dp.lha / dp_main.s < prev    next >
Encoding:
Text File  |  1996-10-18  |  10.4 KB  |  678 lines
  1. ;
  2. ; Double-precision math interface library for ACE.
  3. ;
  4. ; Register a6 is preserved and restored by each routine
  5. ; since Sozobon C uses this as a link register.
  6. ; 
  7. ; Author: David J Benn
  8. ;   Date: 11th February 1996,
  9. ;      16th August 1996,
  10. ;      29th September 1996,
  11. ;      1st,18th October 1996
  12. ;
  13.  
  14. TRUE    equ    -1
  15. FALSE    equ    0
  16.  
  17.     ; Basic functions.
  18.     xref    _LVOIEEEDPAbs
  19.     xref    _LVOIEEEDPAdd
  20.     xref    _LVOIEEEDPCeil
  21.     xref    _LVOIEEEDPCmp
  22.     xref    _LVOIEEEDPDiv
  23.     xref    _LVOIEEEDPFix
  24.     xref    _LVOIEEEDPFloor
  25.     xref    _LVOIEEEDPFlt
  26.     xref    _LVOIEEEDPMul
  27.     xref    _LVOIEEEDPNeg
  28.     xref    _LVOIEEEDPSub
  29.     xref    _LVOIEEEDPTst
  30.  
  31.     ; Transcendental functions.
  32.     xref    _LVOIEEEDPAsin
  33.     xref    _LVOIEEEDPAcos
  34.     xref    _LVOIEEEDPAtan
  35.     xref    _LVOIEEEDPSin
  36.     xref    _LVOIEEEDPCos
  37.     xref    _LVOIEEEDPTan
  38.     xref    _LVOIEEEDPSincos
  39.     xref    _LVOIEEEDPSinh
  40.     xref    _LVOIEEEDPCosh
  41.     xref    _LVOIEEEDPTanh
  42.     xref    _LVOIEEEDPExp
  43.     xref    _LVOIEEEDPFieee
  44.     xref    _LVOIEEEDPLog
  45.     xref    _LVOIEEEDPLog10
  46.     xref    _LVOIEEEDPPow
  47.     xref    _LVOIEEEDPSqrt
  48.     xref    _LVOIEEEDPTieee
  49.  
  50.     ; Miscellaneous external references.
  51.     xref    _AbsExecBase
  52.     xref    _LVOOpenLibrary
  53.     xref    _LVOCloseLibrary
  54.  
  55.     xref    _MathTransBase
  56.     xref    _LVOSPFieee    
  57.     xref    _LVOSPTieee    
  58.  
  59.     ; Exports.
  60.     xdef    _MathIeeeDoubBasBase
  61.     xdef    _MathIeeeDoubTransBase
  62.  
  63.     section    code, dp_code
  64.  
  65. ;************************************
  66.  
  67.     PUBLIC    _dp_open
  68.     ;
  69.     ; Open IEEE DP libraries x 2.
  70.     ;
  71. _dp_open:
  72.     move.l    a6,_a6
  73.  
  74.     movea.l _AbsExecBase,a6        
  75.  
  76.     ; Open basic library.
  77.     movea.l    #_ieeedoubbaslibname,a1
  78.     move.l    #34,d0
  79.     jsr     _LVOOpenLibrary(a6)
  80.     tst.l    d0
  81.     beq.s    _abort_dp_open
  82.     move.l    d0,_MathIeeeDoubBasBase
  83.  
  84.     ; Open transcendental library.
  85.     movea.l    #_ieeedoubtranslibname,a1
  86.     move.l    #34,d0
  87.     jsr     _LVOOpenLibrary(a6)
  88.     tst.l    d0
  89.     beq.s    _abort_dp_open
  90.     move.l    d0,_MathIeeeDoubTransBase
  91.     
  92.     ; Exit with success code (TRUE).
  93.     move.l    #TRUE,d0    ; okay.
  94.     bra.s    _exit_dp_open
  95.  
  96. _abort_dp_open:
  97.     move.l    #FALSE,d0    ; error.
  98.  
  99. _exit_dp_open:
  100.     movea.l    _a6,a6
  101.     
  102.     rts
  103.  
  104. ;************************************
  105.  
  106.     PUBLIC    _dp_close
  107.     ;
  108.     ; Close IEEE DP libraries x 2.
  109.     ;
  110. _dp_close:
  111.     move.l    a6,_a6
  112.  
  113.     tst.l    _MathIeeeDoubBasBase
  114.     beq.s    _exit_dp_close
  115.  
  116.     movea.l _AbsExecBase,a6
  117.  
  118.     move.l  _MathIeeeDoubBasBase,a1
  119.     jsr     _LVOCloseLibrary(a6)
  120.  
  121.     tst.l    _MathIeeeDoubTransBase
  122.     beq.s    _exit_dp_close
  123.  
  124.     move.l  _MathIeeeDoubTransBase,a1
  125.     jsr     _LVOCloseLibrary(a6)
  126.  
  127. _exit_dp_close:
  128.     movea.l    _a6,a6
  129.     
  130.     rts
  131.  
  132. ;************************************
  133.  
  134.     PUBLIC    _dp_longint_to_double
  135.     ;
  136.     ; Convert LONGINT to DOUBLE.
  137.     ;
  138. _dp_longint_to_double:
  139.     move.l    a6,_a6
  140.  
  141.     move.l    8(sp),d0        ; integer value.
  142.  
  143.     movea.l    _MathIeeeDoubBasBase,a6
  144.     jsr    _LVOIEEEDPFlt(a6)
  145.  
  146.     move.l    4(sp),a0        ; address of double result.
  147.     move.l    d0,(a0)+
  148.     move.l    d1,(a0)
  149.  
  150.     movea.l    _a6,a6
  151.     
  152.     rts
  153.  
  154. ;************************************
  155.  
  156.     PUBLIC    _dp_double_to_longint
  157.     ;
  158.     ; Convert DOUBLE to LONGINT.
  159.     ;
  160. _dp_double_to_longint:
  161.     move.l    a6,_a6
  162.  
  163.     move.l    8(sp),a1        ; address of double.
  164.     
  165.     move.l    (a1)+,d0
  166.     move.l    (a1),d1
  167.  
  168.     movea.l    _MathIeeeDoubBasBase,a6
  169.     jsr    _LVOIEEEDPFix(a6)
  170.  
  171.     move.l    4(sp),a0        ; address of integer result.
  172.     move.l    d0,(a0)
  173.  
  174.     movea.l    _a6,a6
  175.  
  176.     rts
  177.  
  178. ;************************************
  179.  
  180.     PUBLIC    _dp_single_to_double
  181.     ;
  182.     ; Convert SINGLE to DOUBLE.
  183.     ;
  184. _dp_single_to_double:
  185.     move.l    a6,_a6
  186.  
  187.     move.l    8(sp),d0        ; single value.
  188.  
  189.     movea.l    _MathTransBase,a6
  190.     jsr    _LVOSPTieee(a6)        ; FFP -> IEEE SP
  191.     movea.l    _MathIeeeDoubTransBase,a6    
  192.     jsr    _LVOIEEEDPFieee(a6)    ; IEEE SP -> IEEE DP
  193.  
  194.     move.l    4(sp),a0        ; address of double result.
  195.     move.l    d0,(a0)+
  196.     move.l    d1,(a0)
  197.  
  198.     movea.l    _a6,a6
  199.     
  200.     rts
  201.  
  202. ;************************************
  203.  
  204.     PUBLIC    _dp_double_to_single
  205.     ;
  206.     ; Convert DOUBLE to SINGLE.
  207.     ;
  208. _dp_double_to_single:
  209.     move.l    a6,_a6
  210.  
  211.     move.l    8(sp),a1        ; address of double.
  212.  
  213.     move.l    (a1)+,d0
  214.     move.l    (a1),d1
  215.  
  216.     movea.l    _MathIeeeDoubTransBase,a6
  217.     jsr    _LVOIEEEDPTieee(a6)    ; IEEE DP -> IEEE SP
  218.     movea.l    _MathTransBase,a6
  219.     jsr    _LVOSPFieee(a6)        ; IEEE SP -> FFP
  220.  
  221.     move.l    4(sp),a0        ; address of single result.
  222.     move.l    d0,(a0)
  223.  
  224.     movea.l    _a6,a6
  225.  
  226.     rts
  227.  
  228. ;************************************
  229.  
  230.     PUBLIC    _dp_add
  231.     ;
  232.     ; Add two DOUBLE values (result = op1 + op2).
  233.     ;
  234. _dp_add:
  235.     move.l    a6,_a6
  236.  
  237.     move.l    12(sp),a2        ; address of operand 2.
  238.     move.l    8(sp),a1        ; address of operand 1.
  239.  
  240.     move.l    (a2)+,d2
  241.     move.l    (a2),d3            
  242.  
  243.     move.l    (a1)+,d0
  244.     move.l    (a1),d1    
  245.  
  246.     movea.l    _MathIeeeDoubBasBase,a6
  247.     jsr    _LVOIEEEDPAdd(a6)
  248.  
  249.     move.l    4(sp),a0        ; address of result.    
  250.     move.l    d0,(a0)+
  251.     move.l    d1,(a0)
  252.  
  253.     movea.l    _a6,a6
  254.  
  255.     rts
  256.  
  257. ;************************************
  258.  
  259.     PUBLIC    _dp_sub
  260.     ;
  261.     ; Subtract two DOUBLE values (result = op1 - op2).
  262.     ;
  263. _dp_sub:
  264.     move.l    a6,_a6
  265.  
  266.     move.l    12(sp),a2        ; address of operand 2.
  267.     move.l    8(sp),a1        ; address of operand 1.
  268.  
  269.     move.l    (a2)+,d2
  270.     move.l    (a2),d3            
  271.  
  272.     move.l    (a1)+,d0
  273.     move.l    (a1),d1    
  274.  
  275.     movea.l    _MathIeeeDoubBasBase,a6
  276.     jsr    _LVOIEEEDPSub(a6)
  277.  
  278.     move.l    4(sp),a0        ; address of result.    
  279.     move.l    d0,(a0)+
  280.     move.l    d1,(a0)
  281.  
  282.     movea.l    _a6,a6
  283.  
  284.     rts
  285.  
  286. ;************************************
  287.  
  288.     PUBLIC    _dp_mul
  289.     ;
  290.     ; Mulitply two DOUBLE values (result = op1 * op2).
  291.     ;
  292. _dp_mul:
  293.     move.l    a6,_a6
  294.  
  295.     move.l    12(sp),a2        ; address of operand 2.
  296.     move.l    8(sp),a1        ; address of operand 1.
  297.  
  298.     move.l    (a2)+,d2
  299.     move.l    (a2),d3            
  300.  
  301.     move.l    (a1)+,d0
  302.     move.l    (a1),d1    
  303.  
  304.     movea.l    _MathIeeeDoubBasBase,a6
  305.     jsr    _LVOIEEEDPMul(a6)
  306.  
  307.     move.l    4(sp),a0        ; address of result.    
  308.     move.l    d0,(a0)+
  309.     move.l    d1,(a0)
  310.  
  311.     movea.l    _a6,a6
  312.  
  313.     rts
  314.  
  315. ;************************************
  316.  
  317.     PUBLIC    _dp_div
  318.     ;
  319.     ; Divide two DOUBLE values (result = op1 / op2).
  320.     ;
  321. _dp_div:
  322.     move.l    a6,_a6
  323.  
  324.     move.l    12(sp),a2        ; address of operand 2.
  325.     move.l    8(sp),a1        ; address of operand 1.
  326.  
  327.     move.l    (a2)+,d2
  328.     move.l    (a2),d3            
  329.  
  330.     move.l    (a1)+,d0
  331.     move.l    (a1),d1    
  332.  
  333.     movea.l    _MathIeeeDoubBasBase,a6
  334.     jsr    _LVOIEEEDPDiv(a6)
  335.  
  336.     move.l    4(sp),a0        ; address of result.    
  337.     move.l    d0,(a0)+
  338.     move.l    d1,(a0)
  339.  
  340.     movea.l    _a6,a6
  341.  
  342.     rts
  343.  
  344. ;************************************
  345.  
  346.     PUBLIC    _dp_cmp
  347.     ;
  348.     ; Compare two DOUBLE values.
  349.     ; Returns an integer:
  350.     ; 
  351.     ;     +1 if op1 > op2
  352.     ;     -1 if op1 < op2
  353.     ;      0 if op1 = op2
  354.     ;
  355. _dp_cmp:
  356.     move.l    a6,_a6
  357.  
  358.     move.l    12(sp),a2        ; address of operand 2.
  359.     move.l    8(sp),a1        ; address of operand 1.
  360.  
  361.     move.l    (a2)+,d2
  362.     move.l    (a2),d3            
  363.  
  364.     move.l    (a1)+,d0
  365.     move.l    (a1),d1    
  366.  
  367.     movea.l    _MathIeeeDoubBasBase,a6
  368.     jsr    _LVOIEEEDPCmp(a6)
  369.  
  370.     move.l    4(sp),a0        ; address of result.    
  371.     move.l    d0,(a0)
  372.  
  373.     movea.l    _a6,a6
  374.  
  375.     rts
  376.  
  377. ;************************************
  378.  
  379.     PUBLIC    _dp_pow
  380.     ;
  381.     ; Raise one DOUBLE to the power of another.
  382.     ; (result = op1 ^ op2).
  383.     ;
  384. _dp_pow:
  385.     move.l    a6,_a6
  386.  
  387.     move.l    12(sp),a2        ; address of operand 2.
  388.     move.l    8(sp),a1        ; address of operand 1.
  389.  
  390.     move.l    (a1)+,d0
  391.     move.l    (a1),d1        
  392.  
  393.     move.l    (a2)+,d2
  394.     move.l    (a2),d3
  395.  
  396.     movea.l    _MathIeeeDoubTransBase,a6
  397.     jsr    _LVOIEEEDPPow(a6)
  398.  
  399.     move.l    4(sp),a0        ; address of result.    
  400.     move.l    d0,(a0)+
  401.     move.l    d1,(a0)
  402.  
  403.     movea.l    _a6,a6
  404.  
  405.     rts
  406.  
  407. ;************************************
  408.  
  409.     PUBLIC    _dp_abs
  410.     ;
  411.     ; Return the absolute value of a DOUBLE.
  412.     ;
  413. _dp_abs:
  414.     move.l    a6,_a6
  415.  
  416.     move.l    8(sp),a1        ; address of operand.
  417.  
  418.     move.l    (a1)+,d0
  419.     move.l    (a1),d1    
  420.  
  421.     movea.l    _MathIeeeDoubBasBase,a6
  422.     jsr    _LVOIEEEDPAbs(a6)
  423.  
  424.     move.l    4(sp),a0        ; address of result.    
  425.     move.l    d0,(a0)+
  426.     move.l    d1,(a0)
  427.  
  428.     movea.l    _a6,a6
  429.  
  430.     rts
  431.  
  432. ;************************************
  433.  
  434.     PUBLIC    _dp_neg
  435.     ;
  436.     ; Return the negated value of a DOUBLE.
  437.     ;
  438. _dp_neg:
  439.     move.l    a6,_a6
  440.  
  441.     move.l    8(sp),a1        ; address of operand.
  442.  
  443.     move.l    (a1)+,d0
  444.     move.l    (a1),d1    
  445.  
  446.     movea.l    _MathIeeeDoubBasBase,a6
  447.     jsr    _LVOIEEEDPNeg(a6)
  448.  
  449.     move.l    4(sp),a0        ; address of result.    
  450.     move.l    d0,(a0)+
  451.     move.l    d1,(a0)
  452.  
  453.     movea.l    _a6,a6
  454.  
  455.     rts
  456.  
  457. ;************************************
  458.  
  459.     PUBLIC    _dp_ceil
  460.     ;
  461.     ; Return the least integer which is >= a DOUBLE.
  462.     ; Result is DOUBLE not LONGINT!
  463.     ;
  464. _dp_ceil:
  465.     move.l    a6,_a6
  466.  
  467.     move.l    8(sp),a1        ; address of operand.
  468.  
  469.     move.l    (a1)+,d0
  470.     move.l    (a1),d1    
  471.  
  472.     movea.l    _MathIeeeDoubBasBase,a6
  473.     jsr    _LVOIEEEDPCeil(a6)
  474.  
  475.     move.l    4(sp),a0        ; address of result.    
  476.     move.l    d0,(a0)+
  477.     move.l    d1,(a0)
  478.  
  479.     movea.l    _a6,a6
  480.  
  481.     rts
  482.  
  483. ;************************************
  484.  
  485.     PUBLIC    _dp_floor
  486.     ;
  487.     ; Return the largest integer which is <= a DOUBLE.
  488.     ; Result is DOUBLE not LONGINT!
  489.     ;
  490. _dp_floor:
  491.     move.l    a6,_a6
  492.  
  493.     move.l    8(sp),a1        ; address of operand.
  494.  
  495.     move.l    (a1)+,d0
  496.     move.l    (a1),d1    
  497.  
  498.     movea.l    _MathIeeeDoubBasBase,a6
  499.     jsr    _LVOIEEEDPFloor(a6)
  500.  
  501.     move.l    4(sp),a0        ; address of result.    
  502.     move.l    d0,(a0)+
  503.     move.l    d1,(a0)
  504.  
  505.     movea.l    _a6,a6
  506.  
  507.     rts
  508.  
  509. ;************************************
  510.  
  511.     PUBLIC    _dp_sin
  512.     ;
  513.     ; Return the sine of a DOUBLE (radians).
  514.     ;
  515. _dp_sin:
  516.     move.l    a6,_a6
  517.  
  518.     move.l    8(sp),a1        ; address of operand.
  519.  
  520.     move.l    (a1)+,d0
  521.     move.l    (a1),d1    
  522.  
  523.     movea.l    _MathIeeeDoubTransBase,a6
  524.     jsr    _LVOIEEEDPSin(a6)
  525.  
  526.     move.l    4(sp),a0        ; address of result.    
  527.     move.l    d0,(a0)+
  528.     move.l    d1,(a0)
  529.  
  530.     movea.l    _a6,a6
  531.  
  532.     rts
  533.  
  534. ;************************************
  535.  
  536.     PUBLIC    _dp_cos
  537.     ;
  538.     ; Return the cosine of a DOUBLE (radians).
  539.     ;
  540. _dp_cos:
  541.     move.l    a6,_a6
  542.  
  543.     move.l    8(sp),a1        ; address of operand.
  544.  
  545.     move.l    (a1)+,d0
  546.     move.l    (a1),d1    
  547.  
  548.     movea.l    _MathIeeeDoubTransBase,a6
  549.     jsr    _LVOIEEEDPCos(a6)
  550.  
  551.     move.l    4(sp),a0        ; address of result.    
  552.     move.l    d0,(a0)+
  553.     move.l    d1,(a0)
  554.  
  555.     movea.l    _a6,a6
  556.  
  557.     rts
  558.  
  559. ;************************************
  560.  
  561.     PUBLIC    _dp_tan
  562.     ;
  563.     ; Return the tangent of a DOUBLE (radians).
  564.     ;
  565. _dp_tan:
  566.     move.l    a6,_a6
  567.  
  568.     move.l    8(sp),a1        ; address of operand.
  569.  
  570.     move.l    (a1)+,d0
  571.     move.l    (a1),d1    
  572.  
  573.     movea.l    _MathIeeeDoubTransBase,a6
  574.     jsr    _LVOIEEEDPTan(a6)
  575.  
  576.     move.l    4(sp),a0        ; address of result.    
  577.     move.l    d0,(a0)+
  578.     move.l    d1,(a0)
  579.  
  580.     movea.l    _a6,a6
  581.  
  582.     rts
  583.  
  584. ;************************************
  585.  
  586.     PUBLIC    _dp_exp
  587.     ;
  588.     ; Return e to the power of a DOUBLE.
  589.     ;
  590. _dp_exp:
  591.     move.l    a6,_a6
  592.  
  593.     move.l    8(sp),a1        ; address of operand.
  594.  
  595.     move.l    (a1)+,d0
  596.     move.l    (a1),d1    
  597.  
  598.     movea.l    _MathIeeeDoubTransBase,a6
  599.     jsr    _LVOIEEEDPExp(a6)
  600.  
  601.     move.l    4(sp),a0        ; address of result.    
  602.     move.l    d0,(a0)+
  603.     move.l    d1,(a0)
  604.  
  605.     movea.l    _a6,a6
  606.  
  607.     rts
  608.  
  609. ;************************************
  610.  
  611.     PUBLIC    _dp_log
  612.     ;
  613.     ; Return the natural logarithm (base e) of a DOUBLE.
  614.     ;
  615. _dp_log:
  616.     move.l    a6,_a6
  617.  
  618.     move.l    8(sp),a1        ; address of operand.
  619.  
  620.     move.l    (a1)+,d0
  621.     move.l    (a1),d1    
  622.  
  623.     movea.l    _MathIeeeDoubTransBase,a6
  624.     jsr    _LVOIEEEDPLog(a6)
  625.  
  626.     move.l    4(sp),a0        ; address of result.    
  627.     move.l    d0,(a0)+
  628.     move.l    d1,(a0)
  629.  
  630.     movea.l    _a6,a6
  631.  
  632.     rts
  633.  
  634. ;************************************
  635.  
  636.     PUBLIC    _dp_sqrt
  637.     ;
  638.     ; Return the square root of a DOUBLE.
  639.     ;
  640. _dp_sqrt:
  641.     move.l    a6,_a6
  642.  
  643.     move.l    8(sp),a1        ; address of operand.
  644.  
  645.     move.l    (a1)+,d0
  646.     move.l    (a1),d1    
  647.  
  648.     movea.l    _MathIeeeDoubTransBase,a6
  649.     jsr    _LVOIEEEDPSqrt(a6)
  650.  
  651.     move.l    4(sp),a0        ; address of result.    
  652.     move.l    d0,(a0)+
  653.     move.l    d1,(a0)
  654.  
  655.     movea.l    _a6,a6
  656.  
  657.     rts
  658.  
  659. ;************************************
  660.         
  661.     section    data, dp_data
  662.  
  663. _ieeedoubbaslibname:    dc.b    'mathieeedoubbas.library',0
  664. _ieeedoubtranslibname:    dc.b    'mathieeedoubtrans.library',0
  665.  
  666. _MathIeeeDoubBasBase:    dc.l     0
  667. _MathIeeeDoubTransBase:    dc.l     0
  668.  
  669. ;************************************
  670.  
  671.     section bss, dp_bss
  672.  
  673. _a6:            ds.l 1
  674.  
  675. ;************************************
  676.  
  677.     END
  678.